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ABSTRACT 

Topological insulators form a novel state of matter that provides new opportunities to create unique quantum phenomena. 
While the materials used so far are based on semiconductors, recent theoretical studies predict that also strongly correlated 
systems can show non-trivial topological properties, thereby allowing even the emergence of surface phenomena that are not 
possible with topological band Insulators. From a practical point of view, it is also expected that strong correlations will reduce 
the disturbing impact of defects or impurities, and at the same increase the Fermi velocities of the topological surface states. 
The challenge is now to discover such correlated materials. Flere, using advanced x-ray spectroscopies In combination with 
band structure calculations, we infer that CeRu4Sn6 is a strongly correlated material with non-trivial topology. 


Introduction 

A metallic surface is one of the necessary characteristics of a 3-dimensional topological insulator. The surface charge carriers 
have a massless Dirac dispersion and are protected from backscattering by time reversal symmetry. So far the materials 
used in the research field of topological insulators are based exclusively on simple semiconductors and only very recently 
it was realized that also strongly correlated systems can be topologically non-trivial. Here the Kondo or intermediate-valent 
insulator SmBg^ came immediately into consideration^"*" and first investigations**"*" show that the robust metallicity of the 
sample surface could be a promising concept for explaining its long-time mysterious low temperature residual conductivity.*^ 
Increasing the interest for correlated topological insulators is the possibility to find phenomena that might be even more 
exotic**"^* than the ones of topological semiconductors. 

The important issue is now to find other strongly correlated materials with non-trivial topology. This is not a trivial 
task since quantitative numerical tools often applied to characterize the topological nature of semiconductors are lacking 
for correlated materials, due to the fact that the reliability and accuracy of mean field solutions of the electronic structure 
in the presence of strong correlations are far from sufficient. It is therefore important that recent theoretical studies have 
drafted effective models for strongly correlated systems that should produce topological non-trivial behavior."’"’*" One of 
the propositions is that rare earth Kondo or mixed-valent insulators should have all the necessary ingredients for realizing 
topological insulating properties:""*" a fully spin-compensated insulating ground state formed by orbitals having opposite 
parity, e.g. 4/ vs. 5d, and inverted in energy in particular sections in the Brillouin zone. It is even specihcally mentioned that 
for Ce systems with cubic"’" and tetragonal symmetry^’ ^ a gap always opens when the ground state is a nodeless 7^ = ±1/2 
Kramers doublet and when the 4/ occupation is sufficiently away from integer. 

Here CeRu 4 Sn 6 comes into play. It is a tetragonal, non-centrosymmetric (/42m),potential Kondo insulator with a 7 = 5/2 
Hund’s rule ground state. A hybridization-induced gap was reported first in 1992 supported further by several experimental 
studies ,24-26 estimated to be of the order of lOOK from single crystal resistivity data .2"’^^ Optical conductivity mea¬ 

surements point towards remnant metallicity and yet the low temperature resistivity is in the milliOhm-cm range, indicating 



clearly the depletion of the charge carriers 23,25,27,28 Magnetization and resistivity are highly anisotropic strongly suggesting 
the J = 5/2 Hund’s rule ground state is split by the tetragonal crystal-electric field (CEF) jxSR measurements did not find 
magnetic order down to 50 mK, thus pointing out the importance of Kondo-type interactions in this material . 

One of the pressing open questions is whether CeRu 4 Sn 6 fulfills the conditions mentioned in the theoretical studies about 
strongly correlated topological insulators To answer this question, we determine here whether the hybridization between 
conduction and 4/ electrons (c/-hybridization) in CeRu 4 Sn 6 is strong enough and whether the CEF ground state has a fa¬ 
vorable symmetry to develop a gap of sizable magnitude. The degree of hybridization has been investigated by performing 
a temperature dependent x-ray absorption spectroscopy (XAS) experiment at the Ce L-edge to obtain the Ce valence and the 
Kondo temperature Tk- The CEF ground state symmetry was studied by measuring the Ce M-edge in an XAS experiment 
and the Ce N-edge in a non-resonant inelastic x-ray scattering (NIXS) experiment. We have also performed band structure 
calculations to understand better our experimental findings and to test the robustness of our conclusions. Details about the 
experimental set-ups and the single crystal preparation and quality can be found in the Section Methods. 

Results 

Cerium 4/ valence and Kondo temperature 

Trivalent cerium (Ce^+) has a 4/ occupancy = 1 (/*). The presence of strong c/-hybridization leads to a partial delocaliza¬ 
tion of the 4/ electrons into the conduction electron sea and the resulting non-integer valent 4/ ground state can be written in 
a local picture as a mixed state |‘Egs) = co |/^) + ci + C 2 |/^L), i-e. with additional contributions of tetravalent (/®) and 
divalent (/^) states.^*’ Here L and L denote the number of ligand holes. For Ce the resulting occupation of the 4/ shell is 
then given by nf = 2c2 + c\ < 1, whereas Cq quantifies the degree of delocalization. The hybridization is described in terms 
of the formation of a Kondo singlet state at ksTx below the energy of the 4/* configuration, which breaks down as thermal 
fluctuations match the same energy. Hence, information about the Kondo temperature Tk is obtained from the temperature 
variation of rif, or equivalently of Cq if C 2 Cj or if C 2 does not change with r. 21-33 

Core level spectroscopy is a powerful tool for determining the occupations of the various 4/" configurations.^*^ This is 
based on the fact that the impact of the core-hole potential depends strongly on the configuration. Moreover, recent XAS 
experiments measuring the Lai recombination process (3d^i2 2 ^ 3 / 2 ) with partial fluorescence yield (PFY-XAS) at the Ce 

Fa-edge {2p —^ 5d) have tremendously improved the accuracy of determining valencies.We have performed such a 
PFY-XAS experiment on CeRu 4 Sn 6 , determined Cq and obtained Tk from its temperature dependence. 

Figure 1(a) shows the PFY-XAS data of CeRu 4 Sn 6 at 5K (open circles). The strongest absorption is due to the /* 
contribution in the ground state, the small but well resolved hump at about 5740 eV is due to the /** contribution, and the small 
shoulder on the low energy tail of the /* peak (just above 5720 eV) is due to some spectral weight. It is not straightforward 
to describe the L 3 -edge spectral line shape since it is determined by the unoccupied 5d density of states in the presence of 
a 2p core hole. Yet, for extracting the /" weights, it is sufficient to fit the data with a set of Gaussian or similar functions 
as a mean to obtain their integrated intensities.We have applied such a fit by describing each spectral weight (/(/**), 
/(/*), and /(/^)) with identical line shapes, consisting of three Gaussian profiles each. Only the total intensities and central 
energy positions were varied. Each profile includes a tanh-type step function (black lines) with the maximum slope at the peak 
maximum to fit the background. The red line in Fig. 1(a) is the result of such a fit. 

For the spectral weights we obtain /(/*’) =0.059 ±0.015, /(/*) =0.931 ±0.015 and /(/^) = 0.011 ±0.015 which corre¬ 
sponds to ~ 0.953 and Cq ~ 0.059 ±0.015 under the assumption !(/") ~ c^. This is an approximation since it neglects 
hybridization effects in the final states.As in core level photoemission,^® 7(/^) overestimates C 2 and /(/**) underestimates Cq 
so that the above values for n f and Cq should be taken as upper and lower limits, respectively. It is therefore safe to conclude 
that the c/-hybridization is very strong in CeRu 4 Sn 6 . 

Figure 1(b) shows the temperature dependence of the PFY-XAS data, normalized to the integrated intensity. The inset 
in Fig. 1(b) shows the /** region on an enlarged scale after smoothing the data, and it is clearly seen that the main changes 
occur within the temperature interval of 150 to 250 K. Applying the same fit by varying only the total intensity of the /* 
and /** contributions and adjusting the height of the background according to the /* absorption line, yields the temperature 
variation Cq{T) [upper inset. Fig. 1(c)]. The dashed line is a broadened step function describing the temperature evolution. 
We determine Tk empirically as the inflection point of c^{T),^^~^^ yielding Lf = (170± 25) K supporting the analysis of the 
temperature dependence of the single crystal resistivity in Ref. 28. 

CEF ground state symmetry 

The local point symmetry at the Ce site of in CeRu 4 Sn 6 gives rise to a tetragonal CEF that splits the J = 5/2 Hund’s 

rule ground state into three Kramers’ doublets; one of the type Fg = 15/2, ±1/2), and two mixed states of type Fj = 
Vl-a2 |5/2, ±3/2) ±a |5/2, t 5/2) and F^ = Vl-a2 |5/2, ±5/2) - a |5/2, t 3/2), in the /-representation We have 
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Figure 1 . (a) PFY-XAS spectra of the Ce L-edge at T = 5 K 
(open circles). The size of the circles indicates the error bars. 
The red line is the total fit, the regions of different gray refer to 
the /(/^) and /(/*') spectral weights, taking into account 
a tanh-type background (black solid lines), (b) Temperature 
evolution of the Ce L 3 -edge spectra and the /(/**) peak region 
on an enlarged scale, (c) Inset with /(/**) ~ Cq as function of 
log(7’). The relative error bars amount to ±0.25%. The black, 
dashed line is a guide to the eye, the red arrow indicates its 
inflection point which is taken as Kondo temperature Tk- 
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Figure 2. Temperature dependence of the 
linearly polarized M-edge XAS spectra of 
CeRu 4 Sn 6 . The arrows mark the spectral 
weight due to the contribution in the ground 
state. 
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Figure 3. (a) Linearly polarized XAS data (dots) at the Ce 
M-edges for r = 4K and simulation with ar6=|5/2,±l /2) 
ground state (lines). Blue corresponds to the polarization 
E II [001] and red to E _L[001]. (b) Experimental (dots) and 
simulated (colored lines) linear dichroism (LD), 
l£±[ooi] ~ l£||[ 00 i]’ of the pure states of the 7 = 5/2 multiplet. 
(c) Same as (b) for the 7 = 7/2 multiplet. The inset in (a) shows 
the spatial distribution of the 4/ electrons for Eg symmetry. 


investigated the CEF ground state wave function with two spectroscopic methods: soft XAS at the Ce M-edge (37 — 4y)40-42 
and hard x-ray NIXS at the Ce N-edge (47 —4/)."*^^ In the soft XAS experiment, the dipole selection rules for linearly 
polarized light give access to the initial state symmetry. Excited CEF states only contribute to the polarization as they get in¬ 
creasingly populated with rising temperature. A NIXS experiment with hard x-rays is truly bulk sensitive, and when performed 
at large momentum transfers \q\ the signal is dominated by higher multipole orders of the scattering cross section.In the 
NIXS experiment the ^-dependence (^-direction with respect to the sample orientation) replaces the polarization dependence 
in XAS and, accordingly, ^-dependent multipole selection rules give access to the ground state symmetry.^^ 

Figure 2 shows the soft XAS data at the Ce M 4 5 edge for the two polarizations E || [001] and E _L[001] for several 
temperatures between 4 and 300 K. At all temperatures there is a strong polarization effect which decreases only slightly 
between 200 and 300 K, thus indicating the first excited CEF state is situated at > 25 meV. The data at 4 K will therefore yield 
ground state information. In addition to the main absorption line, there are smaller humps at slightly higher photon energies 
(see arrows) due to some weight in the initial state. Its spectral weight decreases with rising temperature in agreement with 
the findings from the PFY-XAS data. 

Figure 3(a) shows the linearly polarized low-temperature data (red and blue dots). For the simulation of the spectra a full 
multiplet code is used (see Methods, Analysis). Figure 3(a) shows a simulation assuming a pure |5/2,±l/2) = r 6 ground 
state (red and blue lines), which is the state with rotational symmetry [see inset in Fig. 3(a)]. The simulation reproduces the 
main features of the spectra and their polarization dependence, but the match is by no means perfect. This issue is investigated 
further in the following. 

In Fig. 3(b) the experimental linear dichroism (LD) (/£x[ooi]‘^£||oob black dots) is compared with the simulated LD of the 
three pure 15/2,7^) states (colored lines). The LD of the pure |5/2,±l/2) = r 6 (orange) and also of the |5/2,±3/2) (blue), 
corresponding to a r 7 -symmetry with a = 0, have the correct sign, but the LD of the latter is too small and any |5/2 ,=f 5/2)- 
admixture (green) in a Ty-type ground state would make the discrepancy even larger. Even if we adopt a Fy-type ground state 
of predominantly |5/2,±3/2), the measured LD will be larger than the simulated one. We therefore conclude that although 
the ground state is not a pure r 6 = |5/2,±l / 2 ), it is still mainly of that character. 

A NIXS experiment at the Ce N-edge was performed in order to rule out surface effects being responsible for the unsatis¬ 
factory agreement between the soft XAS data and the simulation. The data in Fig. 4 (red and blue dots) were obtained after 
240 s counting per point and after subtracting a rising, linear background (Compton). The colors are chosen as before, i.e. red 
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Figure 4. Dots: Experimental NIXS intensity at the Ce A ^4 5 
edge after subtracting a linear background. The point size 
reflects the statistical error, the counting time was 240 s per 
point. Lines: Calculated scattering function S(q, CO) for the 
pure states of the J = 5/2 multiplet. Data and simulation are 
shown for q || [ 100 ] (blue) and q || [ 001 ] (red). 


for q II [001] and blue for q || [100]. Again there is a clear dichroic effect. We compare the measured A 4 5 edge with simulations 
for the pure \5/2,J^) states (red and blue lines), which here are shown again (cf. Ref. 43) for reasons of convenience. Details 
about the theoretical simulations containing a full multiplet ansatz can be found in Section Methods, Analysis. 

The strongest dichroism in the data appears at the shoulder structures at about 105 and 113 eV. Qualitatively the simulation 
for the |5/2, ±1 /2) state reproduces the data quite well. In contrast, the |5/2, ±3/2) state suggests a crossing of both directions 
on both shoulders, which is not observed experimentally, and the dichroism of the |5/2,±5/2) state has the wrong sign. We 
therefore conclude also from this experiment that the ground state must be governed by a Te = |5/2, ± 1 /2) symmetry, but at 
the same time is not purely given by the Tg. 


Discussion 

What could be the reason for the deviation between the experiment and the simulation starting from a pure r6 = |5/2,±l/2)? 
Already in 1986 van der Laan et aJ. pointed out that in the presence of strong hybridization higher lying states contribute to 
the M-edge signal In CeRu 4 Sn 6 we hnd that the Kondo temperature (170 K) is not much smaller than the CEF excitations 
(> 25 meV = 290 K) so that a contribution of higher lying CEF states to the ground state is indeed to be expected. Given the 
fact that the LD of these higher lying CEF states is smaller or even of opposite sign, see Fig. 3(b) and (c), it is evident that the 
LD of this composite ground state will be reduced from that of the pure |5/2, ± 1 /2) state. We could also reverse the argumen¬ 
tation: the observation of a reduced LD can be taken as evidence for the presence of strong hybridization. Nevertheless, we can 
conclude that the ground state of CeRu 4 Sn 6 is predominantly of |5/2,±1 /2) symmetry. This is in good agreement with the 
anisotropy of the low-temperature magnetization and susceptibility,^^ thereby supporting also the results of LDAh-DMFT 
calculations at 290 K 

We now will use density functional theory based band structure calculations to obtain a deeper understanding of the 
experimental hndings (details in Section Materials and Methods, Band structure calculations). In Fig. 5 (a) we display our 
result of a full-relativistic non-magnetic calculation employing the local density approximation (LDA) for CeRu 4 Sn 6 . We can 
clearly observe that the band structure has a true band gap in agreement with an earlier study in comprising of unoccupied 
Ce 4/ bands above the Fermi level and occupied Ru 4d bands below the Fermi level. The presence of the gap is important 
since it confirms that the system has the even number of electrons per unit cell necessary to allow for the formation of a fully 
spin-compensated, non-magnetic state. We also notice in the left panel of Fig. 5 (a) that the lowest unoccupied Ce 4/ state is 
surprisingly given by the 7^ = |5/2, ±3/2) (blue dots) and not by the |5/2, ±1/2) (red dots). Yet, looking at the Ce 4/ partial 
density of states depicted in the right panel of Fig. 5 (a) we hnd that it is the |5/2,±1 /2) which has the largest occupation, as 
found in our experiment. We therefore can infer from these calculations that the empty |5/2,±1 /2) has indeed the stronger 
hybridization with the occupied Ru 4d bands so that it gets pushed up higher in energy and, at the same time, acquire a larger 
occupation. 

Next, we will make an attempt to estimate the effect of correlations on the relevant symmetries of the system. We do 
this by introducing artihcially an extra potential for a particular state so that it becomes more occupied, i.e. as to mimic 
the Ce 4/* valence state thereby making the inversion with the Ru 4d. Figure 5 (b) shows the result for moving down the 
|5/2,±1 /2) orbital and Fig. 5 (c) the = |5/2,±3/2). We can observe that in the = |5/2,±3/2) case the system becomes 
metallic with quite a large number of bands crossing the Fermi level For the the |5/2,±l/2) case, we also see that the 
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Figure 5. (a) Full relativistic (LDA+spin-orbit) band structure of CeRu 4 Sn 6 , wherein the Ce 4/ weight is mostly above the 
Fermi level. From the projected DOS, we note the largest contribution to the occupied valence sector is from J = 5/2 and 
= ± 1/2, (b) Applying a constant potential to the = 1/2 state, (c) Applying a constant potential to the = 3/2 state. 


system is metallic, but only with a few crossings and small Fermi surface areas, suggesting that the number of charge carriers 
will be low. Important for the topology is that the band structure is that of a semi-metal, i.e. that there is a warped gap, 
so that a unique assignment can be made about the character of the bands. This aspect did not necessarily follow from the 
conditions set by Dzero and Galitzki,^’ ^ since there could always have been some other non-bonding bands crossing the Fermi 
level, which is not considered in their effective model. Yet, our calculations indicate specifically that such is not the case for 
CeRu 4 Sn 6 . It is also important to notice that the warped gap is substantial, i.e. 20 meV or larger, which corroborates nicely 
with the fact that the 7^ = ± 1 /2 is nodeless and therefore guarantees the presence of hybridization with the Ru 47 bands at all 
k-points of the Brillouin Zone (except at the time reversal inversion points). 

To summarize, we thus have found that CeRu 4 Sn 6 fulfills the conditions mentioned by Dzero and Galitzki^’^ to be a 
strong topological insulator. It has a ground state symmetry given primarily by the Fg = |5/2, ±1 /2) and has a 4/ occupation 
sufficiently far away from integer with a very high Kondo temperature (170K) so that the effective degeneracy is larger than 
two. We have in addition verified that there are no non-bonding bands that otherwise could spoil the formation of a warped 
gap. We therefore conclude that CeRu 4 Sn 6 is a strongly correlated material with non-trivial topology. 

Methods 

Samples 

All experiments were performed on single crystals. Single crystals were grown with the self-flux floating-zone melting method 
using optical heating in a four-mirror furnace. The quality of the crystals was checked with Laue and x-ray diffraction as 
well as scanning electron microscopy and energy dispersive x-ray spectroscopy, which show that the samples are single grain 
and single phase materials. The lattice parameters are in agreement with a = 6.8810A and c = 9.7520A as given in Ref. 22. 
The single crystals were pre-aligned with a Laue camera. A peculiarity of the present tetragonal structure is that y/la Ri c 
within Ri 0.2%. As a consequence the structure appears pseudo-cubic in Laue images and thus the final alignment of the single 
crystals had to be done via magnetization measurements.^^ 
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Figure 6. Dots: isotropic spectrum as constructed from the 
experimental data /^o = 1/3 (/^upoi] + 2/£x[ooi])- Solid line: fit 
with full muliplet routine after adjusting reduction factors (see 
text). 


Set-ups of x-ray techniques 

The PFY-XAS experiment was performed at the GALAXIES inelastic scattering endstation at the SOLEIL synchrotron 
in Erance.^' The synchrotron radiation was monochromatized using a Si(lll) nitrogen-cooled fixed-exit double-crystal 
monochromator {AE/Ef^lAx 10^^), followed by a Pd-coated spherical collimating mirror. The x-rays were then focused to a 
spot size of 30 fim x 90 fim EWHM [vertical x horizontal] at the sample position by a 3:1 toroid Pd-coated mirror. The overall 
energy resolution for the experiment was measured to be 700 meV EWHM with about 1 x lO'^ photons/second on the sample. 
The incident x-ray energy was scanned across the Ce L 3 absorption edge (ft! 5723 eV) while monitoring the intensity of the 
Ce 'La^ de-excitation (ft 4839 eV) using a Ge 331 spherical bent crystal analyzer (Bragg angle 80.8°, 1 m radius of curvature) 
in the Rowland circle geometry which focused the emitted x-rays onto a silicon avalanche photodiode. All measurements 
took place in a backscattering geometry and normalization was performed using a monitor placed just before the sample that 
recorded the scattered radiation from a Kapton foil placed in the beam path. The flight path was hlled with helium gas to 
reduce the air scattering and the samples were cooled in a helium cryostat. 

The XAS spectra were recorded at the synchrotron light source BESSY II using the UE46 PGM-1 undulator beamline. 
The total electron yield (TEY) was measured under ultra high vacuum of the order of 10^**'mbar. The CeRu 4 Sn 6 single 
crystals were cleaved in situ to guarantee clean sample surfaces. The TEY signal was normalized to the incoming photon flux 
Iq as measured at the refocusing mirror. The energy resolution at the cerium M 4 5 edges was set to 0.15 eV. The undulator 
combined with a normal incident measurement geometry allow a change in polarization without changing the probing spot on 
the sample surface. The crystals were oriented such that the two polarizations E || [001] and E _L [001], [001] being the long 
tetragonal axis, could be measured. Several crystals were measured or re-cleaved in order to assure the reproducibility of the 
data. 

The NIXS measurements were performed at the beamline ID20 at ESRE. A single Si(lll) monochromator set the in¬ 
cident energy to 9690 eV and the scattered intensity was analyzed by three columns of three Si(660) analyzers (in total 
nine) at in-plane scattering angles of 20 = 140°, 146°, and 153°, corresponding to an averaged momentum transfer of 
1^1 ft (9.4 ± 0.2) *. This conhguration yields an instrumental energy resolution of EWHM = 1.3 eV and made the Ce N 4 5 

edge measurable beyond the dipole limit.Two crystals with polished (100) and (001) surfaces were mounted in a cryostat 
and measured in specular geometry with q || [ 001 ] and ql. [ 100 ], i.e. parallel and perpendicular to the c axis. 

Analysis 

The soft XAS M-edge data were analyzed with the full multiplet code XTLS 8.3 as described in Ref. 40. The atomic param¬ 
eters are given by Hartree-Eock values and the reduction factors for the Coulomb and spin orbit interaction interaction are 
determined by htting the isotropic spectrum liso = 1/3 (l£||[ooi] +2l£x[ooi]) in Fig- 6 . We obtain reduction factors of 41 % for 
the 4/-4/ and 20% for the 3d-4f Coulomb interaction while the spin orbit values remain unchanged within 1 or 2%. These 
are typical values for Ce intermetallic compounds.Eor the reproduction of the polarized data, the CEE has to be taken 
into account. The calculations are purely ionic. 

The NIXS data were simulated using the full multiplet code Quanty.^^ A Gaussian and a Lorentzian broadening of 
EWHM= 1.3 eV and 0.2 eV, respectively, are assumed to account for the instrumental resolution and life time effects. The 
atomic Hartree-Eock values were adjusted via the peak positions, resulting in reductions of 30 % and 20 % for the 4/-4/ and 
4(7-4/ Coulomb interactions, respectively. The simulation procedure is explained in detail in Ref. 43. 

Band structure calculations 

Density functional theory based calculations were performed using the full-potential non-orthogonal local orbital (FPLO)^^ 
code employing the local density approximation (EDA) and including spin-orbit coupling. The band structure and density of 
states (DOS) are obtained from a carefully converged calculation with 1000 k-points in the irreducible wedge of the Brillouin 
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zone. For easy comparison with literature, we have plotted the band structure along the A:-path as described in Ref. 48. To 
simulate the enhancement in the occupation of the various 7^ states, we have applied a constant potential shift to the on-site 
energies (1 eV) of the relevant orbitals. 
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